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Abstract 

Effective- field-theory methods are used to study the high T limit of QCD. These 
methods unravel the contributions to the free energy of QCD at high temperature 
from the scales T, gT, and g 2 T. The free energy is explicitly computed to order g 5 . 
Implications for the application of perturbative QCD to the quark-gluon plasma are 
also discussed. 



1 Introduction 

Hadronic matter is expected to undergo a phase transition when T ~ T c ~ 200 MeV between 
a low temperature phase in which it is confined in the form of hadrons and a high temperature 
phase in which quarks and gluons are deconfmed. The latter phase, know as Quark- Gluon 
Plasma (QGP), can be studied by perturbative methods since the strong coupling constant 
a s is expected to be small at high temperature. 

A field theory at high temperature may be described by a theory in 3 dimensions resulting 
from decoupling the modes with zero Matsubara frequency 0, |3[]. This dimensionally 
reduced theory can be interpreted as an effective field theory whose parameters are to be 
computed as a perturbative expansion in the coupling constant of the original theory ||]. 

The effective theory approach is now a well-established perturbative method to study 
field theories at finite temperature [|J ^ It is a useful tool for organizing calculations of 
high-order terms in the perturbation expansion and it is especially powerful in dealing with 
non-Abelian gauge theories. It is also an important conceptual tool that explicitly separates 
the different energy scales. Here I present its application to the calculation up to order g 5 
of the free energy []5| and use the separation of scales to analyze the convergence of the 
series [§, |J. 



2 The free Energy of QCD to order g 5 

The partition function of QCD at finite temperature is obtained from the Lagrangian £qcd = 
(1 / '4)Cr" +(77 ^q, with gauge coupling constant g. We consider the effective theory that 
results from integrating out the nonstatic modes. Such a theory, which is called Electrostatic 
QCD (EQCD), is made out of only the static bosonic modes and 3 dimensional. The free 

1 Invited talk presented at Workshop on QCD in Paris, June 1996. 
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energy density is Fqcd — T (fg — \ogZ E Q CD /V), where ^eqcd is the partition function of 
EQCD obtained from the Lagrangian 

£eqcd = \G%G% + i(AA,) 2 + \m\A\ + \\Al + 6C EQCD (1) 

whith effective gauge coupling constant qe- <5£eqcd represents an infinite series of non- 
renormalizable terms. The effects of the fermions are incorporated into the effective parame- 
ters. By power counting, the effective parameters at leading order are: ~ T 3 , m 2 E ~ g 2 T 2 , 
9e ~ 9 2 T, and X E ~ g*T. 

The magnetostatic fields remain massless; therefore, we can go further in separating 
the different scales of QCD at high temperature by integrating out A Q . We obtain an 
effective theory of EQCD which is called magnetostatic QCD (MQCD). The free energy 
density of QCD can be written as F = T(/e + }'m + fa), where fo = —logZyiQCD/V, 
and .Zmqcd is the partition function of MQCD obtained from the Lagrangian £mqcd = 
(1/4)G° G^ + 5£mqcd • The gauge coupling constant is #m and SCmqcd represents an infinite 
series of non-renormalizable terms. 

Again, we can use power counting to identify the order of the leading contribution to 
the MQCD parameters: fu ~ rn E ~ g 3 T 3 and g 2 M ~ g% ~ g 2 T. \ E contributes to the free 
energy only at order g 6 ; so, if we are interested in the free energy at lower order, we can 
ignore Xe- Similarly the non-renormalizable terms of EQCD can also be omitted. 

The leading contribution to fa can be obtained by realizing that the only parameter 
with dimensions involved in MQCD is gu and therefore the leading contribution is of order 
9m ~ 9 6 T 3 . Since we are interested in computing the free energy of QCD up to order g 5 we 
can ignore the contribution from fc; it will be enough to compute Je and fu ■ We conclude 
that the free energy of QCD up to order g 5 is given by F = T[/s(T, g; Ae) + gE', A#)], 

where A^ is a factorization scale that separates the scales T and gT. Therefore, we have to 
determine m 2 E , g E , and fu- Calculations for SU(N) gauge theory with rif fermions and 
results in analytical form are detailed in 0; here, I will just state the results in numerical 
form for SU(3) with 3 fermions. 

The effective mass is computed by matching the electrostatic screening mass for QCD 
and EQCD 

V \ 9 



m E = 0.93 T 



9 2 (H) [ 



4:71 



1 + - 0.26 + 1.43 log- 



(2) 



2-kT ) Ati 

At this order in g 2 , there is no dependence on the factorization scale A^. 

The parameter is determined by calculating the free energy in both full QCD and 
EQCD, and matching the two results. 

f E (A E ) = -14.8 T 3 ' 



1 - 0.90 — 4 

Att 



[18.3 + 6.91 log— - 1.30 log— 
1 5 2vrT & 2ttT/ 
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where g is the coupling constant in the MS renormalization scheme at the scale /i. We have 
used the renormalization group equation of the coupling constant to shift the scale of the 
running coupling constant to an arbitrary renormalization scale \i. 
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Through order g 5 , Jm is proportional to the logarithm of the partition function for EQCD: 
fu = ~ log Z E q C d/V. Now, we have to consider the cotnribution to logZ E Q CD of orders 
g 3 , g 4 , and g 5 which are given by the sum of 1-loop, 2-loop, and 3-loop diagrams. The details 
of this calculation can be found in ||; the final result is 

f M (A E ) = -0.21m| 

The coefficient Jm in (§) can be expanded in powers of g by setting g 2 E = g 2 T and 
by substituting the expression (|2|) for m 2 E . The complete free energy to order g 5 is then 
F = (/e + fhi)T; in agreement with the result obtained independently by Kastening and 
Zhai 0. Note that the dependence on the arbitrary factorization scale Ae cancels between 
J'e and Jm, up to corrections that are higher order in g, leaving a logarithm of T/rriE- 

3 Convergence of Perturbation Theory 

We have calculated the free energy as a perturbation expansion in powers of g to order g 5 . In 
this section, we examine the convergence of the series and it is based on |5], || . We now ask 
how small a s = g 2 / (4-7r) must be in order for the perturbation expansion to be well-behaved 
in the sense that the term of order n is smaller than the term of order n — 1. If the series 
is apparently convergent in this sense, then it can plausibly be used to evaluate the free 
energy. For simplicity, we consider the case rif = 3. If we choose the renormalization scale 
fi = 2nT which is the mass of the lightest nonstatic mode, the correction to the leading 
order result is a multiplicative factor 1 — 0.9a s + 3.3a;?/ 2 + (7.1 + 3.51oga s )«g — 20. 8a^ 2 . 
The a b J 2 term is the largest correction unless a s (27iT) < 0.16 which requires T greater 

than about 15 GeV; at this value of T, the series is 1 - 0.14 + 0.21 - 0.02 - 0.21 H , 

and its convergence is questionable. Now suposse the temperature is a few times above the 
transition temperature (~ 200 MeV); e.g., T ~ 350 MeV then a s (2nT) ~ 0.3 and the series 
has the form 1 — 0.27 + 0.54 + 0.26 — 1.03 + • ■ •. It is clear that the perturbative series is 
not convergent at this temperature. If we consider T ~ 1 TeV, then a s {2itT) ~ 0.07 and 
the series is 1 — 0.063 + 0.061 — 0.011 — 0.027 + • ■ -. However, one has to realize that at such 
high temperatures electroweak processes become important and that we must also take into 
account the effects of heavy quarks. 

We can go further in our analysis provided that we have separated the contributions from 
the scales T and gT. The free energy is the sum of two terms: and ju which contain 
the contributions of order T and gT respectively. The term is given in ([|), when nf = 3 
and = 2tiT, the correction to the leading order result form the series 1 — 0.90a s + (6.47 — 
6.91 log + ■ • ■. The next-to-leading-order correction to Je is independent of jj and 

Ae, and is small compared to the leading-order term provided that a s (2irT) < 1.1 which 
corresponds to T > 67 MeV. The next-to-next-to-leading-order correction can be made small 
by adjusting A E ; it vanishes for A E = 5.1 nT. We conclude that the perturbation series for f E 
is well-behaved if the factorization scale is chosen to be approximately 5ttT. Whether this 
choice is reasonable can only be determined by calculating other EQCD parameters to higher 
order to see if the same choice leads to well-behaved perturbation series. It is interesting to 
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I - (0.54 + 0.721(^7^ ) ( ^ ) -0.70 I 



(4) 
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study the convergence of the other parameter of EQCD that we have computed, m E , which 
is given by (||D and when nf = 3 and /i = 2irT, the correction to the leading order result form 
the series 1 — 0.26a s . For the next-to-leading-order correction to m 2 E to be smaller than the 
leading-order term, we must have a s {2itT) < 3.8, which corresponds to T > 40 MeV. Based 
on these results, we conclude that the perturbation series for the parameters of EQCD are 
well-behaved provided that a s (27rT) < 1, which corresponds to T > 70 MeV. 

We now consider the convergence of the perturbation series (|j) for Jm- The size of the 
next-to- leading-order correction depends on the choice of the factorization scale A E . It is 
small if is chosen to be approximately m E . The next-to-next-to-leading-order correction 
in (|4]) is independent of any arbitrary scales. If rif = 3, it is smaller than the leading order 
term only if a s < 0.17, which corresponds to T > 2 GeV. Thus the perturbation series for 
Jm is well-behaved only for temperatures that are much higher than those required for the 
parameters of EQCD to have well-behaved perturbation series. 

This analysis indicates that the slow convergence of the expansion for F in powers of 
tJoT s can be attributed to the slow convergence of perturbation theory at the scale gT. 

4 Conclusions 

Our explicit calculations were significantly streamlined by using effective-field-theory meth- 
ods to reduce every step of the calculation to one that involves only a single momentum 
scale. They allow us to study the convergence of the perturbation expansion for thermal 
QCD. At the scale T, perturbation corrections can be small only if T > 70 MeV. At the 
scale gT, perturbation corrections can be small only if T > 2 GeV. Thus, in order to achieve 
a given relative accuracy, the temperature must be much larger for perturbation theory at 
the scale gT compared to perturbation theory at the scale T. 

There is a range of temperatures in which perturbation theory at the scale gT has broken 
down, but perturbation theory at the scale T is reasonably accurate. In this case, one can still 
use perturbation theory at the scale T to calculate the parameters in the EQCD Lagrangian. 
However, nonperturbative methods, such as lattice simulations of EQCD, are required to 
calculate the effects of the smaller momentum scales gT and g 2 T. The effective-field-theory 
approach provides a dramatic savings in resources for numerical computation: the effective 
field theory is 3 dimensional and quarks are integrated out of the theory, which reduces it 
to a purely bosonic problem. 
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